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Abstract 

In this paper, we investigate three symmetry breaking effects in strong and radiative decays 
of strange heavy mesons. We study l/mg corrections within the heavy quark effect theory, as 
well as SU{3) and SU{2) symmetry breakings induced by light-quark mass differences and the 
g — TT mixing vertex. These effects are studied in a covariant model. The numerical results show 
that the l/mg corrections of the coupling constants are consistent with a*A qcd • The SU{3) 
symmetry violating effect of the strong coupling constant is obviously larger than that of the 
magnetic coupling constant. The value of the g — n mixing vertex has some changes because of 
the renewed data. As compared with the other theoretical calculations and the experimental data, 
our radiative decay rates are much larger than those of the other theoretical methods, except for 
yPT; however, our branching ratios are close to the experimental data. 
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I. INTRODUCTION 


For excited strange heavy mesons pion and/or photon emissions are the domi¬ 


nant decay modes which determine their lifetimes [l|. Of these decay modes, the radiative 
decay, D* Dsj, and the only kinematically allowed strong decay, D* —>■ -D^vr, which is the 
isospin-violating mode, have been observed, and the branching ratio TJD* Ds'n')/T{D* —)■ 
Dg'f) has been measured by the CLEO and BaBar collaborations . The latter collabo¬ 
ration obtained r(U* —)■ Ds'n')/T{D* —)■ Ds'y) = 0.062±0.005(stat.)±0.006(syst.), which was 
a signihcant improvement over the former one. This precise value provides an ideal occasion 
to test different theoretical estimations for the strong and electromagnetic interactions of 
strange heavy mesons. 

In 1989, it was realized that, in low-energy situations where the typical gluon momenta 
are small compared with the heavy quark mass (mg), quantum chromodynamics (C^D) 
dynamics becomes independent of the flavor (mass) and spin of the heavy quark 
These new spin and flavor symmetries combine to form an SU{2NJ) symmetry, called 
heavy quark symmetry (HQS), which is not manifest in the original QCD Lagrangian. HQS 
allows us to factorize the complicated light quark and gluon dynamics from that of the 
heavy one, and thus provides a clearer physical picture in the study of heavy quark physics. 
Beyond the symmetry limit, a heavy quark effective theory (HQET) can be developed by 
systematically expanding the QCD LaCTangian in powers of I/ttiq, with which HQS breaking 
effects can be studied order by order jo-lsj. Although the development of HQET from QCD 
has simplihed the analysis of heavy hadron physics, many properties of hadrons, for example, 
their decay constants and axial coupling constants, are still not calculable directly from QCD. 


To study these quantities, one unavoidably has to use phenomenological models 


the structures of hadrons. These include the constituem 
MIT bag model [ll, the lattice QCD calculations 
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quark model (CQM) 


;o describe 
1^, Q , the 


14| . QCD sum rules [1^, and the 


light-front quark model (LFQM) lfil-ll8|. In spite of the fact that the CQM and the MIT bag 
models have been widely used, the results calculated from these two models are trustworthy 
only for processes involving small momentum transfers. The LFQM is a relativistic quark 
model with simple boost kinematics which allows us to describe physical processes with large 
momentum transfers. However, this model is not a fully Lorentz covariant [^, and this 
defect limits its usefulness to matrix elements with space-like momentum transfers < 0) 
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o„.y. Moreover, the LFQM is not eapable of handi.ng the so-called Z-diagrams M- In 
Ref. [19|, a covariant light-front model of heavy mesons has been suggested. However, 
the approach taken there is not systematic, and light-quark currents are not considered. 
To overcome the drawbacks mentioned above, a covariant held theoretical model has been 
proposed for the heavy meson bound state problem 


23l |. This model is fully covariant 
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and satishes HQS; at the same time, it retains the simplicity of the quark model picture. 
This theory allows us to formulate theoretical calculations in terms of the standard Feynman 
diagrams. Therefore, the lack of Z-diagrams in the ordinary LFQM is no longer a problem. 
Combining this model with HQET, we can systematically study various l/uig corrections 
to heavy meson properties in the framework of perturbative held theory. 

In the other extreme, due to the relatively small light-quark masses the 

light-quark sector of the QCD Lagrangian obeys an approximate SU{3)l x SU{3)r chiral 


symmetry 


24| . Due to the spontaneous breaking of the chiral symmetry, there exist eight 


pseudoscalar bosons (called Goldstone bosons, which include three tt’s, four iF’s, and one 
p), whose dynamics obeys the SU{3)l x SU{3)r chiral symmetry. If we want to study the 
low-energy interactions of heavy hadrons and Goldstone bosons, we need to build an effective 


theory that obeys both chiral and heavy quark symmetries. This was done in references 


25 


29| , where chiral symmetry and HQS were synthesized in a single effective chiral Lagrangian 


which described the strong interactions between heavy hadrons and Goldstone bosons. The 
theory has since been extended to incorporate electromagnetic interactions as well 281-182 1. 
In principle, the effective chiral Lagrangian provides an ideal framework in which to study 
the strong decay mode. However, symmetry considerations alone, in general, do not lead 
to quantitative predictions, unless further assumptions are made to extract the values of 
the various coupling constants appearing in the Lagrangian. Furthermore, the framework of 
an effective chiral Lagrangian does not allow for a systematic discussion of HQS violating 
l/mg effects, which is important for a thorough understanding of heavy quark physics. In 
fact, in the heavy-light {Qs) system, there are three different types of symmetry breaking 
mechanisms: (1) HQS breaking from l/mq corrections, ( 2 ) SU{3) symmetry breaking due 
to strange quark mass (m^ ^ mu^), and (3) SU{2) symmetry breaking due to the up-down 
quark mass difference {niu 7 ^ m^). The purpose of this paper is to systematically study 
these symmetry breaking effects in a covariant model for the strong and radiative decays of 
strange heavy mesons. 
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The paper is organized as follows. In Sec. II, we briefly review the covariant model, 
which is based on HQET. Some heavy meson properties in the heavy qnark limit and l/mg 
corrections are considered. The nnmerical calcnlations and discnssions are expressed in Sec. 
III. In Sec. IV, we make some conclnding remarks. 

II. FORMALISM 

The covariant model starts from HQET in the heavy qnark limit (mg oo) and describes 
a heavy meson as a composite particle, consisting of a rednced heavy qnark conpled with a 
brown mnck of light degrees of freedom. It is formnlated in an effective Lagrangian approach, 
so that it is fnlly covariant, and we can nse Feynman diagrammatic techniqnes to evalnate 
varions processes. 


A. Covariant model 


Using the l/rnq expansion to the heavy qnark QCD Lagrangian g, [Tj], the QCD La¬ 
grangian for heavy and light qnarks pins glnons can be written as L = Lq + L^q , where 


Lq = h^iv ■ Dhy + q - rUg) q - -F^''Fa^u, 


LrriQ — 


E 

n=l 


2m, 


Q 


hyi p±{—iv ■D)"' ^i]/)±hv, 


( 2 , 1 ) 

( 2 . 2 ) 


— v^v ■ D is orthogonal to the heavy qnark velocity, Lq is responsible for binding a 
heavy qnark and a light qnark in the heavy qnark limit, and Lmn contains l/mq corrections 


to Lq. The effective Lagrangian we 
of psendoscalar heavy mesons reads 


rave constrncted to describe the low-energy dynamics 
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23|: 


Lefi ^ r + d + h.c. 


(2.3) 


where represent the composite psendoscalar heavy meson fields which appear only as 
external states, 

A = lim mju — (2-4) 

rriQ^oo 

is their residnal mass in the heavy qnark limit. 


q^ = GF{—iv ■ d)q 


(2.5) 
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FIG. 1: Feynman rules in the heavy quark limit. 


represents collectively the degrees of freedom in a heavy meson, where F is a form factor 
whose presence is expected for an effective interaction resnlting from the non-pertnrbative 
QCD dynamics, and G is the normalization constant given by 


G-^ = i 


d^p 


F‘^{v ■ p) 


V ■ p + niq 


(A — n • pY^p^ — mV) 


( 2 . 6 ) 


At this point, we note that F{v-p) is analogons to the meson wave fnnction in the LFQM, and 
G is the corresponding normalization constant. To explicitly evalnate G and other physical 
qnantities, we need to specify the structnre fnnction F{v ■ p), which is nnfortnnately not 
calcnlable from hrst principles. Nevertheless, from the constraints that F does not depend 
on the heavy qnark residnal momentum and it forbids on-shell dissociation of the heavy 
meson into Qq, a plausible form for F is: 


F{v ■ p) = ip{v ■ p){A — V ■ p), 


(2.7) 


where the function ip{v ■ p) does not have a pole at n • p = A. 

Within this framework, hadronic matrix elements are calculated via standard Feynman 
diagrams where heavy mesons always appear as external legs. The Feynman rules for this 
effective theory are shown in Figure [H Figure [U^a) specihes the meson-Q-g vertex with 
Tm = ^ 75 (~ /); for M is the pseudoscalar (vector) meson. All the other Feynman rules 
are the same as in QCD and HQET. Thus, if the power of \p\ in the wave function p is 
less than — |, this model will work well. After building a covariant framework to describe 
heavy meson structures, we go on to evaluate some of the basic heavy meson properties. 
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FIG. 2; Feynman diagram for heavy meson decay constant. 

These include the decay constant, the l/mq corrections of the heavy meson mass, and the 
axial-vector and electromagnetic coupling constants of strange heavy mesons. 


B. Decay constants and I/ttiq corrections of the heavy meson mass 


Consider the heavy meson decay constants dehned by: 

{^\qi^l5hv\P{v)) = ifpv^, 

mrK\V{e)) = /veC 


The Feynman diagram to be evaluated is illustrated in Figure |2J Using the Feynman rules 
in Figure [H the matrix element is evaluated as: 


(0|V-,r,ft„|M(.)) = f A GF(v.p)(v.p^ra,) ^ 


= /MTr 


(27r)^ (p2 — rn'^ y . p 


— '^)Tm 


( 2 . 8 ) 


where AI, = 3 is the number of colors, while y/Nc arises from the color wave function of 
meson, (rf-|- gg + bb)/\/Nc, and Tm = il 5 {— ^ for a pseudoscalar (vector) heavy meson; 
the corresponding weak current vertex is F^ = 7 ^ 75 (y^^). Here, as mentioned in the last 
subsection, the meson held is represented by the form factor F. Thus, the decay constant 
in the heavy quark limit is given by: 


fM = 2y/3iG 


(Pp ip{v ■ p){v ■ p + niq) 
(27r)^ {p^ — rn?q) 


(2.9) 


We hnd that this decay constant is the same for pseudoscalar and vector heavy mesons, 
which is in accord with the prediction of HQS. Jm is related to the usual dehnition of decay 
constant fM by fM = fiM/y/mM- 
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FIG. 3: (a), (b), and (c) are Feynman rules for and (d) is the light quark coupling 

to a gluon. 


Next, we consider the l/mg corrections of the heavy meson mass. The Lagrangian in 
Eq. (I2.2|l can be expanded as: 

LmQ = 0i + 02 + 0 > ( 2 - 10 ) 

where 0^ = ^K K, O 2 = ^ K G,, K, and = T^Gf = 

is the gluon field strength tensor. Oi is the gauge invariant extension of the kinetic energy 
arising from the off-shell residual motion of the heavy quark, and O 2 describes the color 
magnetic interaction of the heavy quark spin with the gluon field. It is clear that both 
Oi and O 2 break the flavor symmetry, while O 2 breaks the spin symmetry. The Feynman 
rules for these HQS breaking interactions are given in Figure O With the l/mg corrections 
included, the heavy meson masses can be expressed as: 

rriM = rng -|- A — --(Ai -|- dM^2), (2-11) 

2mQ 

where (Im = 3(—1) for the pseudoscalar (vector) meson, Ai comes from Oi, and A 2 comes 
from 02 - Ai receives two different contributions, which are a kinetic energy piece and a 
one-gluon exchange piece, thus, Ai = A) ^ A) A The relevant Feynman diagrams are shown 
in Figure m Using Feynman rules in Figures [1] and [3l we can readily write down the various 
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FIG. 4: Feynman diagrams for IjmQ corrections to meson mass. 


contributions: 


= iG^ 


A<»' = -C,G\ 
<^mA2 = —glCfC"^ 




(27r)"^ (p2 — rri^ 

2 2 f d^P d^P’ 


2(p^ — V ■ p‘^){v ■ p + mq), 

(/3^(n • p')(p{v ■ p) 


{2tiY{2tiY (p'2 — rriY){p^ — rng){p — p') 
d^p d^p' pYv ■ p')(p{v ■ p) 


{2'kY{‘^t^Y {p'"^ ~ Tkn?){p‘^ — m'i){p — p 


_ 7"i 

/'\2 OMi 

r2 

I AA 


/A 2 


where G/ = | is a color factor and TjP are defined by: 

Tm = 2{{p ■ p' + p'‘^ — V ■ pv ■ p' — V ■ p'Y{'mq + V ■ p) 

+ {p ■ p' + p^ — V ■ pv ■ p' — V ■ pYi'mg + V ■ p')}, 

4 

Tm = -g1m{(p - P ■ p' + V ■ pti ■ p' - V ■ P )(m, + V ■ p) 

-{p-p' — p — V ■ pv ■ p' + V - p‘^)(mq + V ■ p)}. 


( 2 . 12 ) 

(2.13) 

(2.14) 


(2.16) 

(2.16) 


As expected, and Xj are the same for both pseudoscalar and vector mesons. The 
hyperhne mass splitting is obtained : 

2 A 2 


= my — mp = 


miQ 


(2.17) 


C. Strong coupling constant 


First, we study the zero order of strong coupling constants. An effective Lagrangian of 
pseudoscalar jP) and vector {V) mesons and their couplings to the Goldstone bosons is 


constructed as 


M- 


Cvp = V>^P^ - MlPP^ + ifMniP - V^A>^PY - 






(2.18) 
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FIG. 5: Feynman diagram of /q. 

where + V^)P\ and V^{A^) is the (axial) vector held: 


A = - ?94*). 

^ is dehned as .^ = ^ M is a 3 x 3 matrix for the octet of Goldstone bosons: 


(2.19) 

( 2 . 20 ) 


M = 


i!_ _|_ JL 

v/2 ^/6 

7r“ 

K- 


TT’’ 


K+ 


.Zh + JL /vO 

V2 ^ vF ^ 


A^' 


-\/i^ 


( 2 . 21 ) 


and /jr is the pion decay constant. Throngh the partial conservation of axial-vector current 
(PCAC), a soft pion amplitude can be related to a matrix element of the axial-vector current 
-4)) = as: 


{BTT-iqm = ^-^{B\Al\A), 

In 


( 2 . 22 ) 


From the chiral Lagrangian, we obtain: 




(2.23) 


On the other hand, the matrix element on the right hand side of Eq. (I2.22p can be evaluated 
in the covariant model. The Feynman diagram to be evaluated is illustrated in Figure [5l 
and the relevant Feynman rules are illustrated in Figure [H The result is: 


(M'(n)|^,T“ 7 ^ 75 ^/>,|M(n)) = Tr 


P ( 1 + '^)-p 

7At75l M' --1 M 


yt Y 

^ ^ AT ^ 


where the x are SU{3) wave functions of the heavy mesons and 


^ 3^^ j {271 


d'^p , , ,, 2 /T 3m‘^ + + 2{v ■ p)"^ + 6mqV ■ p 

- imv-p)\ {A-v-p) - 2 - 


(p2 _ 


(2.24) 


(2.25) 
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FIG. 6: llmq corrections to strong / and magnetic d coupling constants. Fq, stands for an external 
current. Other notations are defined in Figure [3l 


For V —)■ Pn (Fm = — /, Fm' = i'Jb), we compare Eq. (\2.24h with Eq. fl2.18p and where 
the subscript 0 denotes zeroth order in I/ttiq. 

Next, we shall calculate the hrst order l/rriQ corrections of strong coupling constants. The 
relevant matrix elements are collectively illustrated in Figure |6l For V —?■ Pn and Fq = Fg], 
we can readily write down the matrix elements and calculate the traces for Figures [6] (a) as: 

ip{v ■ p)(p{v ■ p') 


“ 2mQ J (27r)'^(27r)^ — m‘^){p'‘^ — m‘^y{p' — pY 


n\ 


(a) 


where 


= \{p'-p + p^ — v-p'v-p—{v- pY) 


{mg + V ■ p'Y + - V ■ p'Y 

o 


+ (p'^ + p' ■ p — (n • p') — V ■ p'v ■ p) - {rUg + V ■ p) {2mg + v ■ p') 


(2.26) 


(2.27) 


Similarly, we can evaluate jW'a for Figure [6] (b), and it turns out that M.a^ = M.^a ■ Then 
a comparison with the chiral Lagrangian result shows: 

glp\v ■ p')p{v ■ p) 


= -2 


CfG^ r d^pdy _ 

niQ J (27r)"‘(27r)"^ (p^ — mPjip'"^ — w?Yip' ~ P) 


'fp (“) 

-/Cl . 


(2.28) 


The above calculation can be repeated for Fq = Fq 2 - We hnd that 5f2 is given by: 


^/2 = 


CfG^ f d^pd^p' 


qIpYv ■ p')f{v ■ p) 


mq J {27rYi‘2TrY {p^ — m?g){p'‘^ — m'^YiP' ~ P) 


np {“) 

2'^2 ) 


(2.29) 
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where 


= 


rl / / f\2 I / 

p — p ■ p — [V ■ p ) + V ■ pv ■ p 


p' ■ p — p‘^ — V ■ p V ■ p + {v ■ pY 


{2m? + 2mgV ■ p + v ■ p' v ■ p — p' ■ p) 


(m^ + 2mqV ■ p + {v ■ p) ) 


(2.30) 


Figure [6] (c) corresponds to the contribution from the heavy quark kinetic energy. For 
V —>■ Ptt, the matrix element can be simplihed as: 


= 5/f ^Tr 


7„75(- /) 


Ai\ 


(2.31) 


where 


^ = 2^^^ 


{p^ — V ■ p^) 


{m + v- pY + -(p^ — n • p^) . (2.32) 
3 


j (27r)"^ (p2 — + ie)^ 2mQ 

Therefore, we obtain the strong coupling constant, including the l/mq corrections, as: 


/ = fo + Sfi^^ + SfY^ + 6f2. 


(2.33) 


D. Magnetic coupling constant 


We now consider the coupling constant which governs the d eca.v V —)■ Py. The relevant 
lowest-order chiral and gauge-invariant Lagrangian is given by 32l |: 


C'yp = Mh X 


-d{eQ? + ?Qe) 


P^^PI’ + h.c.. 


(2.34) 


where 


Q = 


i 0 0 

0 f 0 
0 0 f 


(2.35) 


is the light-quark charge. In the mg oo limit, the Feynman diagram to be calculated is 
similar to Figure El except that the axial-vector current ^4“ is replaced by the light-quark 
electromagnetic current = eeqYq^fj.'Yq- The result is: 


{M'{v)\Yq{ieeq^Y'^q\M{v)) = VeqTi 


Jr ^-p 

—Fm 


(2.36) 
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where 


V = 2ieG^ j ■p)WK-v-p) > (2-37) 

and g = p' — p —)■ 0. For V —)■ P 7 {Tm = — /, F^' = ilb), we compare Eq. fl2.36p with Eq. 
( 12 .34^ and obtain do = V12. 

Next, we calculate l/mg corrections to the magnetic coupling d corresponding to E —?■ 
F 7 . The relevant Feynman diagrams are shown in Figure [ 6 ] with Fq, = ieeq^a- For V P 7 , 

{P'y{q,e)\V{e)) = iCq 2d , (2.38) 


which comes from the effective chiral Lagrangian Eq. fl2.34p . the calculated procedures 

are similar to those of the strong coupling constants. Here, we only show the results for 
p _pis) 


^^{g) ^ f d^p dV (pl(n • p')^{v • p)5f°^ 


Smg J (27r)® (p2 — m2)(p'2 — mg)2(p'— p)2’ 


(2.39) 


where 


= - 2 < 


p^^ + p^ ■ p — (n • pO — V ■ p'v ■ p 


p' ■ p + p^ — V ■ p' V ■ p — {v ■ pY 


(n • p + rriq) 


{vp' + mq)\, (2.40) 


and for Fq = F^, 


hdo = 


25 'sG^ f d'^P d'^p' pl(n • p')<p(n • p)S: 


(a) 


Sut-q j (27r)® (p2 — mH(p'2 — m2)2(p'— p)^’ 


(2.41) 


where 


d“> = i 


/2 ! t l\2 I ! 

P — p ■ p — [V ■ p ) + V ■ pv ■ p 


(n • p + TUq) 


p' ■ p — p^ — V ■ p' V ■ p + {v ■ pY 


(n • p' + rriq) 


For Figure E] (c), we obtain: 


^ _.^^2 f d^p ifjv -pYiv-pP mq){p^ -v-p^) 


(27r 


2mQ(jP‘ — mYY 


(2.42) 


(2.43) 
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In radiative decay, there is an additional l/mg correction which comes from the magnetic 
moment of the heavy quark. The matrix element of this process is: 


= —G^ [ ^ \ (n • p + mg)ieQe^^«;3g^n°6^,(2.44) 


mt 


2mQ''^^" tuq J (27r)4(p2 

if Tm' = * 75 , Tm = — From the normalization condition given in Eq. (I2.6p . we obtain: 

e 


SdQ = 


2m, 


Q 


(2.45) 


Including the above results, we can write: 


d = do + ^ + 6d[^^ + Sd2, 

dq = d -\ — —dq. 
eq 


(2.46) 


III. NUMERICAL RESULTS AND DISCUSSION 


For obtaining numerical results, we shall further assume the form of (p(n • p): (i) ip{v ■ p) 
is an analytic function apart from isolated singularities in the complex plane, and (ii) it 
vanishes as |n • p| oo. These two conditions allow us to evaluate the (or p~—) 
integrations in Eq. fl2.6p by Cauchy’s Theorem. Thus, we take: 

1 


^n{v-p) = , . . 

[y ■ p^ u — le) 

which was used in a previous work 


{n = integer). 


(3.1) 


33| . There are some parameters {ms,mQ,u,as) in this 


covariant model, and we follow the strategy described below to £x them. In a quark model, 
flavor SU{3) symmetry is broken because the strange quark mass (m^) is quite different 
from the up or down quark mass {mu,d)- However, the size of the difference. 


6mq = m 


s 


(3.2) 


is not accurately known. For current quark masses, the value of 5mq was quoted as Smq{p = 
1 GeV) ~ 190 MeV 341 and Smq^p = 2 GeV) ~ 90 MeV in the different renormalization 
scales. On the other hand, for constituent quarks in a relativistic quark model, one typically 


gets 


35|: 


6mq ~ 140 ~ 200 MeV. 


(3.3) 
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Because Srug is an important parameter in our calculations, using the variant values in the 
above range will leading to the quite different results. This will slash our predictive ability. 
Here we use the constraint that is independent of the heavy quark mass to obtain the 
value of SrUg. In other words, Srriq is no longer a free parameter. The processes are as follows: 
we first quote the value mu,d = 0.245 GeV from the previous work and try the initial value 
with 6mq = 140 MeV. Subsequently, we take the charm quark mass and the quark-gluon 


coupling to be the same as that for the non-strange charm meson 33|], and choose an u to 
calculate and A 2 from Eqs. fl2.12p — l|2.14p . Using Eq. fl2.17p . the value of lj can 


be adjusted to fit the hyperfine mass splitting 




AMd*Ds = 143.8 ± 0.4 MeV. 


(3.4) 


After fixing u, we take the bottom quark mass and the quark-gluon coupling to be the same 
as that for the non-strange bottom meson 33| to estimate the other hyperfine mass splitting, 
AMbjBs- In addition, using Eq. (12.lip , we can determine two values of for both the 
charm and bottom sectors. Because A^ is independent of the heavy quark mass, the above 
processes are repeated by fine-tuning the value of 6mq until the two values of A^ are the 
same. Finally, the decay constant in the heavy quark limit, /m^, can also be evaluated in 
terms of Eq. 02.91) . These results are listed in Tables [H and Ull 


n 

6mq (GeV) 

mq (GeV) 


w (GeV) 

A 2 (GeV2) 

Ai (GeV^) 

As (GeV) 

/m. (GeV3/2) 

8 

0.225 

1.73 

0.400 

1.19 

0.124 

-0.210 

0.290 

0.507 

10 

0.219 

1.72 

0.392 

1.78 

0.124 

-0.232 

0.285 

0.508 

12 

0.215 

1.72 

0.387 

2.38 

0.124 

-0.244 

0.281 

0.508 


TABLE I: Us-meson parameters for 


n 

6mq (GeV) 

mq (GeV) 

Cts 

a; (GeV) 

A 2 (GeV2) 

Ai (GeV^) 

AMb*b, (MeV) 

As (GeV) 

8 

0.225 

5.09 

0.381 

1.19 

0.118 

-0.230 

46.3 

0.290 

10 

0.219 

5.09 

0.373 

1.78 

0.118 

-0.251 

46.3 

0.285 

12 

0.215 

5.09 

0.368 

2.38 

0.118 

-0.264 

46.3 

0.281 


TABLE II: B^-meson parameters for 
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First of all, we see that the choice of ipn {n = 8,10,12) makes very little difference. The 
value of Srriq = 215 ~ 225 MeV is close to the typical light-quark mass, Eq. fl3.3|) . used in 


a relativistic formalism 


351]. The hyperhne mass splitting AMb*Bs is consistent with the 


average data: = 46.1 ± 1.5 MeV. However, the value of ag in Bg meson seems to 

be rather larger than the one which is determined by the perturbative evolution equation 
(at the one-loop level in the MS scheme): 




1 -F ag{Mz)/3o\n[{ml°^''/Mz)'^]/{4:7r] 


~ 0 . 22 , 


(3.5) 


where mf'" = 4.89 GeV, Mz = 91.19 GeV, /^o = 11 - fAT = 11 - | for Nf = 4, and 


i\f - 

tf- 


ag{Mz) = 0.119 from experimental hts. The reason is that [7[ if the gluons which coupling 
to the heavy quarks are hard (i.e., the virtual momenta is of order of the heavy quark 
mass), they can resolve the nonlocality of the propagator of the small component helds 
Hy = i p±h^/{iv ■ D + 2mQ — ie). Their effects are not taken into account in the naive 
operator product expansion which was used in the derivation of the effective Lagrangian 
in Eq. fl2.10p . Thus, HQET provides an appropriate description only at scales fi <C rrih, 
and the relevant ag in the 5-quark mesons will larger than ag{ml°^‘^). As to the reduced 
mass, we compare = 0.281 ~ 0.290 GeV with that of the non-strange heavy meson, 
A = 0.202 ~ 0.210 GeV 33], and hnd that the residual mass difference is only about 80 
MeV, in contrast to Srriq = 215 ~ 225 MeV. This can be understood as follows. Due to its 
heavier mass, the strange quark is more tightly bound than an up or down quark; thus, part 
of the mass difference 6mq is compensated for by a larger binding energy of the (Qs)-system. 
We can then obtain the predicted meson decay constant by using = fMs/\lM ms 
and the ratio fMs/fM- 


/m, ^219 MeV, 


IMs 

Jm 


= 1.13 ±0.05, 


(3.6) 


where the value fM — fB = 194 ± 9 MeV (an average of t 
is chosen. For comparison, the QGD sum rules results of 

Ib. 


le results 


36, 


37 ] in lattice QGD) 


fB, = 242+11 MeV, 


/s 


_ 11 7 + 0.03 

~ ' - 0.045 


(3.7) 


and 


39 


4Q| 


fp 

Jb, = 225.6 ± 18.3 ± 3 MeV, ^ = 1.184 ± 0.023 ± 0.007, 

JB 


(3.8) 
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(3.9) 


and the lattice QCD calculation results ^ of 


fB, = 224(5) MeV, ^ = 1.205(7), 

JB 

are shown here. In Tables I and II, the kinetic and chromomagnetic expectation values, Ai 
and A 2 , are the heavy-strange meson parameters which were dehned in some papers as /i^ 
and respectively. The relations between them are 42|: 


Ai — —A2 — /^g/3- 


(3.10) 


Although these heavy-strange meson parameters are not found in the other theoretical cal¬ 
culations, here we show the relevant parameters of B meson [n = 10) which obtained from 
the previous work 33|: Ai = —0.162 GeV^,A 2 = 0.117 GeV^. Gompared with the recent 


result which comes from the inclusive decays with = 4.553 GeV (Eq. fl3.10l) is 
used): Ai = —0.465 GeV^,A 2 = 0.111 GeV^, and we hnd that A 2 is consistent with ours 
because of a bound from the B hyperhne splitting, but otherwise Ai is rather different from 
ours. 

Next, the l/rriQ corrections to fs and dg are listed in Table HTTl In order to do a comparison 



DtDs 

BtBg 

n 

fsO dfsl Sfs2 fs 

dgo ddgx ^ds2 dg 

dgQ ^ds\ ^ds2 dg 

8 

-3.29 0.292 -0.190 -3.18 

0.529 -0.0382 0.0307 0.522 

0.529 -0.0146 0.00988 0.525 

10 

-3.10 0.303 -0.184 -2.98 

0.505 -0.0399 0.0299 0.495 

0.505 -0.0150 0.00962 0.500 

12 

-2.99 0.308 -0.181 -2.86 

0.491 -0.0408 0.0294 0.479 

0.491 -0.0153 0.00948 0.485 


TABLE III: l/mg corrections to / and d. 

for the l/rriQ effects, we introduce an effective gluon mass of mg ~ Aqcd — 300 MeV. For 
D*Ds mesons, 

= 5.97 %, = 5.88 %, 

are consistent with the rough estimate of 

^ (6 ~ 7) %, (3.11) 


16 






















v 

Z2(V) 



(a) 
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FIG. 7; Renormalizations of (a) the strong coupling constant / and (b) the FP 7 coupling constant 
d in chiral perturbation theory. 


for as = 0.4. For BlBg mesons, since mb/rric — 3, consequently, the HQS violating effect and 
the rough estimation are both smaller than those for the D*Ds system by approximately a 
factor of 3. Additionally, as compared with that of the non-strange charm mesons 33|, we 
see that the SU{3) breaking is severe for fs, but less so for dg (see Table HVll . The reason for 


mu,d{GfN) 

/ 

d(GeV)-i 

(5mg(GeV) 


d,(GeV)-i 

0.245 

-1.13 

0.361 

0.225 

-3.18 

0.522 


TABLE IV: SU{3) symmetry breakings to / and d for n = 8 . 


this can be traced back to the fact that the strong coupling constant is sensitively dependent 
on As, but the magnetic coupling constant is sensitively dependent on both A^ and rris in 
this model. The details are as follows. From Eqs. f|2.25p and fl2.37|) . both the strong and 
magnetic coupling constants are dependent on A*. As the power of \p\ in the wave function 
(p must be smaller than — | (see Eq. fl2.6p ). the strong and magnetic coupling constants 
satisfy a simple relation, dso = which is similar to the result in Appendix A of Ref. 

. Combined with the estimation that the total l/mq correction is about (3 ~ 4)% for 


the D*Ds mesons, we obtain an approximate equation, dg — Therefore, the SU{3) 

breaking of dg has been reduced by the factor nig in the denominator. 

The study of SU{3) breaking in chiral perturbation theory follows a different route, in 
which SU{3) symmetry is assumed at the tree level and symmetry breaking effects are 
induced via meson loops (see ^ for details). Thus, from Figure [7] (a) we have: 


/ = 


^ZAV)ZAP)ZM) ,0 


ZiiVPct)) 


r 


(3.12) 
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where Z\ and Z 2 are, respectively, the wave function and vertex renormalization constants, 
(j) denotes a Goldstone boson, and /° is the unrenormalized coupling constant. The Zs have 
all been evaluated in 29|]. Putting in the numbers in fl3.12p . we obtain: 


/ = 1.33, fs = 1.47, 


for /° = 0.52, which hts to the experimental data for non-strange mesons 29j. Thus, we see 
that in chiral perturbation theory, SU{3) breaking in the strong coupling constant is not 
large, with: 


fs-f 


f 


0.10 


(3.13) 


chiral 


This is very different from what we found in the covariant model. As for the radiative decay 
constants in chiral perturbation theory, we have (see Figure [7] (b)): 


d = 




(3.14) 


Z,(VPt) 

where do is the unrenormalized transition magnetic moment, and = 0.394 GeV“^ is 
obtained from htting to the branching ratios of D* —>■ D'y 29|]. Putting the numbers in 
(I3.14p . we obtain d = 0.436 GeV“^, ds = 0.575 GeV“^ and 


d. — d 


d 


= 0.319. 


chiral 


The latter one is close to that of our model: 

d. — d 


d 


= 0.446. 


Finally, we consider the decay widths T{D* — )■ DgTi^), T{D* —>■ Ds'j), T{B* —)■ and 

the ratio: 

t{d: ^ oy) 


r., = 


T{Dt ^ ZI.7) ’ 

which is known to be = 0.062 ± 0.008 experimentally [l|. Note that the decay mode 


(3.15) 


D* ^Ds + tt^ 


(3.16) 


violates isos 
46| 


order 


45 


pin or SU{2) symmetry, and it must proceed via rj — n mixing in the leading 
, as depicted in Figured (a), where = (7r°|i7em|h) is the rj-Ti mixing vertex. 
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(a) 


77" 

/ 


/ 


'eff 

(b) 


Dc 


FIG. 8: D* Ds + TT^ via y^vr-mixing mechanism. 


Figure [HI (a) can be replaced by Figured (b) with an effective D*Ds7i coupling constant: 


-^eff — 


H, 


r}7r 


(3.17) 


y/s m2 - m2 

The strength of the rj-n mixing interaction, ifryTr, can be calculated in various models 
Here we chose to utilize the experimental rates of Svr®. From the data of Particle Data 
Group |i|, we can use F(p —)■ all) = 1.31 ± 0.05 keV and Br(p —)■ 37 r°) = (32.68 ± 0.23)% to 
obtain: 


F(p ^ 37 r°) = 0.428 ± 0.019 keV. 


(3.18) 


fl 


As to the amplitude of // —)■ dvr", a fit of the data in Ref. [l[ shows to be essentially 

constant over phase space: = Mq(1 + 2az), where z is the square of the relative 

distance to the center of the Dalitz plot and a = —0.0315 ± 0.0015. Then the three-body 
phase space integral for constant amplitude was estimated in Ref. [48|: 


F (?7 —)■ 37 r°) ~ 0.827 keV. 

Combining Eqs. fl3.18p and flO.lOp . we obtain the constant amplitude: 


(3.19) 


|M^^ 3 ^o| = 0.719 ±0.033. 


(3.20) 


On the other hand, from the current-algebra PCAC 49 


summing the three cyclic permutations of Figured 50|, 


the total amplitude of n —)■ 37r° is 


M 


-77^377^ 


(7r°|j7em|h) 
m2 - ml 


mt 


2-^r]^27r'^) - M^\r] ^ r]l27r^) 

J 77 


0 (7r°|i7em|h') ml-ml 


(7r0|i7em|h) ml, - ml 


(3.?1) 
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(b) 
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(c) 


FIG. 9: One of three cyclic amplitudes in 77 —)• Svr^ 


where the hrst, second, and third terms of Eq. (I3.2ip correspond to Hg in Figs. [9](a), (b), and 
(c), respectively. The former is the Weinberg tttt tttt ^ strong amplitnde extrapolated 
to the T] mass shell consistent with four-momentum conservation, and the latter is the strong 
amplitude rl^^2'K^) of the 77 ^'^ pole which extrapolated from the 77 ^'^ mass to the vr 

mass. Summing the cyclic permutations and considering the 7777 ' mixing, the authors of Ref. 


5l| obtained: —)■ 77 ,r 2 vr°) = 6 cos ^0 where 0 is the mixing angle: 


I 77 ) = cos 01 7777 ) — sin 01 ss), 
I 77 ') = sin 0 |f 7 r 7 )cos 0 |ss), 


(3.22) 


with I 7777 ) = (I 7777 ) -|- \dd))/\/2, and the third term of E g. fl3.2ip is negligible. The mixing 


angle can be determined by a theory of particle mixing 


m- 


tan 0 = 


{ml, - 2m\ + ml){ml - ml) 


(3.23) 


(2777- 2777^ - 7772 ) ( 777 ^, - 
from which one obtains: 0 = 41.5°. Combining Eqs. fl3.20p and fl3.2ip with the pion decay 
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n 

D*g DgTT^ 

Dg —>■ Dgy 

rg 

Dg —> Dg'y 

8 

0.277«;°i(7.29l°:«^%) 

3.53(92.7 =F 0.7%) 

0.07861°;°°™ 

0.407 

10 

0.243;0;°24(7.48l°J?%) 

3.00(92.5 =F 0.7%) 

0.0809l°;r77 

0.371 

12 

0.224;0;0i (7.601°;™%) 

2.72(92.4 =F 0.7%) 

U.UOZO_Q QQ79 

0.351 


TABLE V: Predicted decay rates (in unit of keV), the branching ratios (in parentheses) and for 

in¬ 
constant /jr = 130.41 ± 0.20 MeV and the relevant meson masses [l|, we obtain: 


= -6400 ± 310 MeVl 


(3.24) 


This value is different from that of Ref. 


^|: = —5900 ± 600 MeV^, the chief reason 


being that the new experimental data in Eq. fl3.18p is used. Then, using the /* and dg 
obtained in Table IIIIl we can predict the ratio and the relevant decay widths in Table El 
For the different in, fhe deviations of the decay rate are about 5 ~ 6 times larger than those 
of Tg. The main reason is that the decay rate T{D* —)■ Dg'y), for example, can be simplihed 
as: 


T{Dl Dg^) cx 


— dg + 


nir 


(3.25) 


where the minus sign comes from the charge of the s quark. On the other hand, the ratio 
Tg can be simplihed as: 

2 


Tg oc 


dgUig 


-dg + ^ 
* rUc 


(3.26) 


where dg ~ is again applied. In other words, the deviations of Vg for the different 

in are suppressed strongly because dg is approximately proportional to fg. Finally, we list 
the predicted decay rates and branch ratios within this work (n = 8) and some theoretical 
models in Table IVIl For comparison, the experimental data are also included. We can End 
that, hrst, our branching ratios and Vg are close to the experimental data. Second, there 
are many theoretical calculations of T{D* Dg'y) in the literature. However, except for 
the result of yPT, our result is much larger than that of the other theoretical methods. In 
fact, the estimations of T{D* —)■ Dg'y) are also quite different among the other theoretical 
groups. A similar situation exists for T{B* — Dg'y). Third, in contrast to T{D* Dg'y), 
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Reaction 

D*g DgTT^ 

Dg Dg^ 

Dg —>• total 

rs(xl0 2) 

Bg Bgj 

Exp.UJ 

(5.9 ±0.7 %) 

(94.2 ± 0.7 %) 

< 1900 

6.2 ±0.8 


This work 

0.277lg;[]|(7.3±0.7%) 

3.53(92.7 =F 0.7%) 

3.56 ±0.03 

7.86l°;?9 

0.407 

XPT [29]t 


4.5 




LFQM [^1 


0.18 ±0.01 



0.068 ±0.017 

RQM 

0.0197 ± 0.0070 (input) 

n “^91 +0.009 

U.OZi_o 008 

0.341 


0.136 ±0.012 

QCDSR 


0.59 ±0.15 




NJLM 


0.09 



0.10 

LQCD 

0.0040 (input) 

0.066 ± 0.026 

0.070 ± 0.028 



MIT 





0.0510 

NRQM 


0.21 




NRQM 


0.40 



0.18 


TABLE VI: Predicted decay rates (in units of keV) and branch ratios (in parentheses) of some 
models. For comparison, the experimental branching ratios are given in the first row. (yPT: chiral 
perturbation, LFQM: light-front quark model, RQM: relativistic quark model, QCDSR: QCD sum 
rules, NJLM: Nambu-Jona-Lasinio model, LQCD: lattice QCD, MIT: MIT bag model, NRQM: 
non-relativistic quark model.) f: The value for g = 0.52, /3 = 2.6 GeV“^, and rric = 1.6 GeV. p. 
The values correspond to a linear model. (): The value for = 0.55. p. The values correspond to 
(a) model. 


there are few cornputations for T{D* — DsTt) violating the isospin symmetry. The relevant 
results of Refs. [5^, come from taking the experimental branching ratio as the input. 
Then, for T{D* Dgir) as well as r^, we need to make further comparisons by means of 
more experiments and theoretical calculations. It is worth mentioning that, in our model, 
although the deviations of are smaller than those of the decay rate for the different (pn 
(because of dg (x fs] see above), this does not mean that we can obtain almost the same no 
matter what, for example, the value of T{D* —)■ Dg'y) is. The average value of T{D* — )■ Dg^) 
for the other theoretical calculations (except yPT) is about one-tenth of ours. From Eqs. 
fl3.25p and fl3.26p . we hnd that if one adjusts the dg to reduce the value of T{D* —)■ Dgj) to 
one-tenth, the value of Vg will be enhanced to about 3 ~ 4 times that of the experimental 
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data. In other words, the fact that our is close to the experimental data gives us conhdence 
in our results and the validity of our covariant framework. 


IV. CONCLUSIONS 


Based on HQET, we have discussed the strong and radiative coupling constants of strange 
heavy mesons in l/mq corrections and SU{3) symmetry breakings. These effects were 
studied using a fully covariant model. The covariant model starts from HQET in the heavy 
quark limit and describes a heavy meson as a composite particle, consisting of a reduced 
heavy quark coupled with a brown muck of light degrees of freedom. It is formulated 
in an effective Lagrangian approach, so that it is fully covariant, and we used Feynman 
diagrammatic techniques to evaluate the various processes. 

The parameters of this model, and u, were chosen to £t the data of the hyperhne 
mass splitting, AMd^Ds^ and because the residual mass is independent of the heavy 
quark mass. Then the other hyperhne mass splitting AMb*Bsj A^ and the decay constant 
in the HQ limit can be calculated. Our AMb*Bs was consistent with the data. The residual 
mass difference between A^ and A was only about 80 MeV, and obviously smaller than that 
between the s and u, d quarks. This is understood as follows. Due to its heavier mass, 
the strange quark is more tightly bound than an up or down quark; thus, part of the mass 
difference between the s and u, d quarks is compensated for by a larger binding energy of 
the (Qs)-system. The SU{3) symmetry breaking effect fj uJ J m is close to the fBs/fs of the 


QCD sum rules 


38 


40l | and the lattice QCD calculation ^ results 


The l/mg corrections of A and ds are consistent with asAQCD/mQ for both the D*Ds 
and the B*Bs systems. In the charmed meson sector, the HQS violating effects are larger 
by approximately a factor of 3 because of mh/rric — 3. The SU{3) symmetry violating 
percentage of fg, at about 180%, is obviously larger than that of the dg at 45%. The reason 
for this was that, due to an approximate equation, dg ~ SU{3) breaking of the dg 

was reduced by the factor rUg in the denominator. For comparison, we estimated the SU (3) 


29| and obtained about 10% 


symmetry violating percentages in chiral perturbation theory 
and 32% for fg and dg, respectively. 

In order to calculate the decay rate of D* —)■ DgTi, which violates isospin or SU{2) 
symmetry, we used the new data of r(77 37r°) to estimate the r] — tt mixing vertex: 
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Hjjn = —6400 ± 310 MeV^. Combining the coupling constants fs, dg, and Hn-K, we studied 
r(Zl* —)■ DgTT^), T{D* —)■ Dg'y), T{B* —)■ Bg'y), and the ratio r^. For the different (fn, the 
deviations of T{D* —)■ DgH^) and T{D* —)■ Dg'j) were about 5 ~ 6 times larger than those 
of Tg. The main reason was dg ~ other words, the deviations of Vg are 

suppressed strongly because dg is approximately proportional to /«. 

Finally, we compared our results with the experimental data and the other theoretical 
calculations in Table VI. Our branching ratios of F(Z1* —)■ DgTi^), T{D* Dg'-)), and were 
close to the experimental data. However, our predictive decay widths of T{D* —)■ Dgn^), 
D* Dg'-), and B* —)■ Bg'-) were much larger than those of the other theoretical groups 
except for yPT. Because another computational Vg is not found in the literature, we tried to 
reduce our T{D* —)■ Dg'-)) to the average value of the other theoretical calculations by directly 
adjusting dg, and hnd that our Vg is enhanced to about 3 ~ 4 times of the experimental data. 
In other words, the fact that our Vg is close to the experimental data gives us conhdence not 
only in the validity of our covariant framework, but also in our predictions about the decay 
widths of D* —)■ DgTi, D* —)■ Dg'-), and B* —)■ Bg'-). Then more experiments about the above 
decay widths are needed. 
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